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ABSTRACT 

We study Giant Molecular Cloud (CMC) environments surrounding 10 Infrared Dark Clouds 
(IRDCs), using ^^CO(l — 0) emission from the Galactic Ring Survey. We measure physical properties 
of these IRDCs/GMCs on a range of scales extending to radii, R, of 30 pc. By comparing different 
methods for defining cloud boundaries and for deriving mass surface densities and velocity dispersions, 
we settle on a preferred “CE,r,G” method of “Connected Extraction” in position-velocity space plus 
Gaussian fitting to opacity-corrected line profiles for velocity dispersion and mass estimation. We 
examine how cloud definition affects measurements of the magnitude and direction of line-of-sight 
velocity gradients and velocity dispersions, including associated dependencies on size scale. CE,r,G- 
defined GMCs show velocity dispersion versus size relations a oc which are consistent with the 
large-scale gradients being caused by turbulence. However, IRDCs have velocity dispersions that 
are moderately enhanced above those predicted by this scaling relation. We examine the dynamical 
state of the clouds finding mean virial parameters dvir — I.O for GMCs and 1.6 for IRDCs, broadly 
consistent with models of magnetized virialized pressure-confined polytropic clouds, but potentially 
indicating that IRDCs have more disturbed kinematics. CE,r,G-defined clouds exhibit a tight corre¬ 
lation of oc with n 0.7 for GMCs and 1.3 for IRDCs (c.f., a value of 0.5 expected for a 

population of virialized clouds). We conclude that while GMCs show evidence for virialization over 
a range of scales, IRDCs may be moderately super virial. Alternatively, IRDCs could be virialized 
but have systematically different ^^CO gas phase abundances, i.e., due to freeze-out, affecting mass 
estimations. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - Stars: formation 


I. INTRODUCTION 

Galactic star formation resides mostly within Giant Molecular C louds (GMCs), conventionally defined to have masses 
> 10^ Mq and observed to extend up to several xlO^ Mq (e.g.. Blitz 1993 Williams et al.|[2QQQ McKee & Ostriker 
2007). With typical mass surface densities of ^ 100 Mq pc“^, GMCs have mean radial sizes of ^ 6 — 100 pc, assuming 


simple circular symmetry. However, GMCs are highly irregular and hierarchical structures. Their dense clumps can 
spawn stellar clusters and associations, creating the bulk of the Galactic field star population. T he efficiency and rate 


of star formation from these clumps is relatively low, i.e., a few percent per local free-fall time (Zuckerman & Evans 
I974j Krumholz fc Tan|2007 ). This appears to be mostly because much of the C MC material is sta ble with respect to 


gravitational collapse, especially material below a threshold Ay ^ 10 mag (e.g., Lada et ^ 2010). Higher total star 
formation effici encies, ^ 10 ~ 50%, appear to be possible in the star-forming clumps that form at least moderately 
bound clusters ( Lada fc Lada||2003 ). 


Different theoretical models of the processes that create star-forming clumps within GMCs or prevent overdensities in 
the bulk of the cloud a re actively debated. These proces ses include the regulation of star formatic)n and stabilization of 
gas b y magnetic fields ( McKee|1989[ Mouschovias)2001 ) or turbulence ( [Krumholz fc McKee|2005 Padoan fc Nordlund 


201 Ij), and/or the initiation of star form ation by d iscrete triggering events, such as convergin g atomic flows (e.g 


Heltsch et al.||2006 ), cloud collisions (e.g., Tan||200Q ) or stellar feedback (e.g., Samal et al.||20i4 ) 


Egan et al.||1998 [Carey et al.[[2000| [Rathborne et al.[[2006| [Butler fc Tan[[2009| Peretto fc Euller|2009| |Battersby et al. 

2010[). Thus their study may help us understand the processes that initiate star formation in GMCs. There h^ 

many investigations of the internal properties of IRDCs, including their temperatures (e.g., [Pillai et al.[[2006 

we been 

Peretto 

et al. 

[20101 Ragan et al.[[2011 

Chira et al. 2013[), mass surface density structure (e.g., Butler fc Tan||2009 

Peretto 

& I’u 

ller||2009 Ragan et al.||20l 

T1 Butler fc Tan [2012 Kainulainen fc Tan[[2013l Butler et al.|2014[), kinematics (e.g.. 


bntani et al.||2006||Hernandez et ^[2011 ), chemistry ( jbanhueza et al.||26l3p . See jTan et al. T fuiir Tor 


CO depletion 
a review. 

However, there have been fewer studies connecting IRDCs to their larger-scale environments, such as the morphology. 
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kinematics and dynamics of their parent clouds. Theories involv ing production of dense g;as in shocks have been 


suppo rted by detection of large-scale SiO emission along IRDCs (Jimenez-Serra et ah 2010 Nguyen-Lu’o’ng et al 


2013|). However, these studies focus only on a few individual clouds and are still conhned to a few-parsec scales in and 
around the filamentary molecular clouds. 

Here we study the ^^CO(l -0)-emitting gas in a nd around 10 well-studied IRDCs, utilizing data from the BU-FCRAO 
Galactic Ring Survey (GRS; Jackson et al.|2006 ). We consider a range of scales out to 30 pc projected radius, expected 
to encompass the potential GIVIC enviromnent of the IRDC. While our study c onnects to scales typical of other large- 
sample GMC studies (e.g., |Heyer et al.||2009| |Roman-Duval et al.||2009| |2010|), by focusing on just 10 regions we are 
able to investigate their kinematic properties in much greater detail. Other studies d one on such a range of scales, 
from clump to GMGs. have been perform ed on nearby GMGs, such as Orion A (e.g., iShimaji ri et al.||2011|), Taurus 
(e.g.. Goldsmith et al.|2008), and Perseus (Ridge et al.|2i 


2006 


iFoster et al.|2009| |Kirk et al.|2010p. However, t 
1 star-forming clumps that is found in IRDGs. 


ese local 


GMGs do not seem to give rise to the more extreme range of 

The main focus of this paper is to “bridge the gap” between IRDC and GMC studies. The questions we aim to 
address include: Are IRDCs typically found within GMGs? Are IRDCs found within specific locations with respect to 
GMCs? Are IRDCs and their surrounding GMCs virialized, and does their degree of virialization vary as a function 
of cloud physical scale? 

The IRDC/GMC sample is presented in §2. Methods for defining cloud boundaries and estimating masses and 
kinemat ic pr operties are described in §3. Results are presented in 5T including derived physical properties of the 
clouds (^.1), the location of I RDCs within their respective GMCs (^4.2), and the kinematic and dynamical analysis 
of the IRDC and GMCs (p.3p. We conclude in §5. 


2. THE IRDC SAMPLE AND MOLECULAR LINE DATA 

We utilize the ^^CO(J = 1-^0) data of the GRS survey ( Jackson et al.[2QQ6 ), which covers 18° < I < 56°, 5±1° and 
—5 < i;lsr < 135 kms“^ (for I < 40°) and —5 < i^lsr < 85 kms“^ (for I > 40°). The GRS had a spatial resolution 
of 46", with 22" sampling, and a spectral resolution of 0.212 kms“^. The typical rms noise is cr(T^) = 0.13 K (or 

^ 0-26 K with a main beam efficiency of 0.48). _ _ __ 

Our selected IRDGs ar e the 10 clouds from [Butler fc Tan| ([2009[ hereafter BT 09) (see also [Butler & Tan||2012t 
[Kainulainen & Tan[[2013[). This sample, a subset of that ol|Ratnborne et al.[ ([2006[), was chosen while consi dering the 
8/im (IRAC band 4) images from the Spitzer Galactic Legacy IVlid-Plane Survey Extraordinaire (GLIMPSE; [Benjamin 


et al.|20Q3). These IRDGs were selected for being relatively nearby, massive, dark (i.e., relatively high contrast against 


the surrounding diffuse emission), and surrounded by relatively smooth dif fuse emission within the 8/im GLIMPSE 
images. Gharacteristic sizes and boundaries for each IRDG were taken from Simon et al. (2006a), where ellipses were 
fitted based on extinction in MSX images. Although these fitted ellipses are not necessarily accurate of IRDC shapes, 
they provide a convenient measure of the approximate cloud structure. The catalog coordinates and sizes for the 10 
IRDCs (including “single component” sub-classifications, labelled “-s”, see below) are listed in Table 


TABLE I 

The IRDC Sample 
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18.822 

-0.285 

14.1 

3.90 

74.0 

18.850 

-0.276 

64.3 

18.862 

-0.276 

64.5 

4.80 

8.90 

B 

19.271 

0.074 

8.40 

1.80 

88.0 

19.287 

0.075 

26.5 

19.287 

0.075 

26.7 

2.40 

7.11 

C 

28.373 

0.076 

12.0 

9.30 

78.0 

28.370 

0.062 

79.0 

28.364 

0.062 

79.2 

5.00 

7.86 

D 

28.531 

-0.251 

20.4 

5.10 

60.0 

28.567 

-0.233 

78.6 

28.567 

-0.239 

80.7 

5.70 

7.00 
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28.555 

-0.251 

87.9 

28.549 

-0.257 

87.7 



E 

28.677 

0.132 

14.6 

4.10 
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7.00 
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34.428 
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57.5 

3.70 

6.48 
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58.0 

34.428 

0.240 

58.0 
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2.00 
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-0.559 

43.9 
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-0.559 

44.1 

2.90 
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34.770 

-0.559 

44.4 

34.764 

-0.559 

44.4 
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35.395 

-0.336 

20.6 

6.40 

59.0 

35.422 

-0.319 

47.5 

35.416 

-0.319 

48.0 

2.90 

6.61 

H-s 






35.391 

-0.337 

43.9 

35.391 

-0.337 

44.4 



I 

38.952 

-0.475 

7.50 

3.00 

64.0 

38.942 

-0.473 

43.5 

38.942 

-0.473 

43.2 

2.70 

7.24 

J 

53.116 

0.054 

1.60 

1.30 

50.0 

53.114 

0.056 

21.4 

53.120 

0.056 

21.6 

1.80 

7.05 


IRDC coordinates and elliptical sizes ad opted from [Simon et aL] ([2006a[) 
.. .... rtr-r: M 


°IRDC kinematic distances adopted from 


Rathborne et al.| (120061 


115 - 

‘^The ^^CO excitation temperatures adopted from Roman-Duval et aF |( [2010[ > 


Kinematic distances were adopted from Rathborne et al. (2006) and Simon et al. (2006b), where the IRDC central 


velocity was matched morphologically between the mid-infrared (MIR) extinction estim ated from MSX a nd the ^^CO 
emission from the GRS. These distances were estimated assuming the rotation curve of Clemens (1985). We assume 
uncertainties of 20 %, given the size of streaming motions of ^10 — 20 kms“^ of clouds with trigonometric parallax 
measurements (e.g.. Brunt haler et al.[[2009 Reid et aT][2014 ). 
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Fig. 1. — The IRDC/GMC sample: each row presents a single cloud (A through E) showin g data o ut to a radius of 30 pc. Left to right: 
GLIMPSE 8 /im image with the IRDG boundary marked with a pink ellipse (see Figs. [ 3 ] and )Al|A4| for intensity scales); GO(1-0) GRS 
integrated intensity map (over the range vq ± 15 kms“^) in units of K kms“^; Full projections of the cloud in position-velocity space: 
simple-extraction (SE) ^^GO(l-O) vi^j. vs I map, SE ^^GO(l-O) vi^j. vs b map; connected extraction (GE) ^^GO(l-O) vi^j. vs I map; GE 
GO (1-0) visr vs b map. The IRDG is marked with a white box of width equal to the extent of the elliptical boundary and height equal 
to the FWHM of the IRDG co-added spectra. The solid black line shows the mass-weighted linear velocity gradient (dvo/ds) across each 
cloud at R = 30 pc. 


3. METHODS: DERIVATIONS OF MOLEGULAR GLOUD PROPERTIES 
3.1. Definition of Cloud Boundaries 

We first determined the cloud center-of-mass in position-velocity space. To do this, we co-added the spectra 

within an elliptical boundary defined by twice the r^aj aiid ^min values, using the full GRS velocity range. Then, a 
velocity interval of 30 km s“^ was considered, centered on the peak of the emission profile. The cloud’s center-of-mass 
was evaluated for both optically thin and opacity corrected gas, (lo^bo^vo) and (/o,r 5 ^o,r 5 '^o,r): respectively. We then 
use this as the reference point around which to search for GMC-scale gas, out to a projected radius of 30 pc and over a 
new velocity range of vq ± 15 km s“^. This range is wide enough to sample the velocities of GMG gas, even those that 
are potentially undergoing collis ional interactions at velocities set by Galactic shear at the tidal radii of the clouds 
(Tan 2000 Tasker & Tan )2009|). The GLIMPSE S/am and GRS ^^GO integrated intensity maps centered on these 


locations are presented in Figures and 
We examined emission within the vq ± 15 km/s velocity range and within apertures of varying radii from the 

cloud center-of-mass coordinate. To explore how boundary definitions affect estimated physical properties, we used 
two different methods to select line emission associated with the cloud. First, “Simple Extraction (SE)” selected 
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Fig. 2.— The same as Figure but for Clouds F through J. 


all the emission within radii i? = 5,10, 20, and 30 pc and vq ± 15 km/s. Second, “Connected Extraction (CE)” 

defined a cloud as a connected structure in /-6-u-space with all cloud voxels required to be above a given threshold 
intensity: we defined a given voxel as “molecular cloud gas” if its CO(1-0) line intensity Tb^v ^ 1-35 K (i.e., the 

GRS ScTrms noise level). The CE search was also limited to within a 30 pc radius and ±15 kms“^ of the IRDC 
center-of-mass. Position-velocity maps for each cloud, defined by both SE and CE, are shown in Eigures[^and|^ 
Clouds defined by CE do not have simple radial sizes. Therefore, we estimated three circular boundaries, centered 
on the extracted cloud’s center-of-mass: (1) mass-weighted radius, Rm, defined as mean projected radial distance of 
cloud mass from the center of mass; (2) areal radius, Ra^ defined by the total projected area A = = NpAp^ where 

Np is the total number of pixels subtended by the cloud and Ap is the area of one image pixel; (3) half-mass radius, 
Ri/ 2 , defined as the radius from cloud center that contains half of the total mass. 


^ 1/2 5 

For 


or the fflPCs, we al s o selec ted “cloud” material via SE and CE using the spatial coordinates and elliptical bound¬ 


aries from Simon et al. (2006a), along with the vq ± 15 kms ^ velocity intervals 


3.2. Column Densities and Masses from ^^CO emission 

We estimated the ^^CO column density of each molecular cloud voxel, dA/'isco from their J = 1 ^ 0 line emission 
assuming a partition function with a thermal distribution described by an excitation temperature Tex via: 


dWsco STrQrot gi 


dv 


Qul- exp{-hiy/[kTeA ’ 


( 1 ) 
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where Qrot is the partition function, A = 6.294 x is the Einstein coefficient, Aq = 0.27204 cm, gi = 1 and 

gu = ^ are the statistical weights of the lower and upper levels and is the optical depth of the line at frequency 
i.e., at velocity v. Each GRS voxel has a velocity width of dv = 0.212 kms“^. Eor linear molecules, the partition 
function is Qrot = + l)exp(—Ej/ZcTex) with Ej = J(J + l)hB where J is the rotational quantum number 

and B = 5.5101 x 10^^ s“^ is the rotational constant. Eor ^^CO(l-O) we have Ej/k = 5.289 K. 

Many studies of the physical properties of GMCs and IRDCs have accounted for line optical depth when estimating 


2011 


their physic al pr operties (e.g., Heyer et al.||2009 Roman-Dnval et ^|2010 Hernandez fc Tan]|2011 Hernandez et al 


2012). In Hernandez & Tan 


(| 2011 |), we showed that optical depth correction factors can increase the ^"^GO 
column density by a factor of ^ 2 m the densest, sub-parsec scale clumps of IRDGs. However, for the more diffuse 
GMGs the optical depth correction factors are expected to be smaller. To gauge the importance of this effect, we carry 
out column density estimates for both the optically thin assumption and accounting for opacity corrections. 

The optical depth is evaluated via 


Tb,u = ^[/(Tex) 


/(Tbg)][l-e- 


( 2 ) 


where Tb,u is the brightness temperature at frequency f{T) = [ex.p{hu/[kT]) — 1] and Tbg = 2.725 K is the 
background temperature. Tb^u is derived from the antenna temperature, Ta^ via Ta = ^/ciump^B,zy, where g is the 
main beam efficiency {g = 0.48 for the GRS) and /dump is the beam dilution factor of the ^^GO emitting gas, which 
we assume to be unity due to the large scale extent of GMGs. Smaller scale structures are undoubtedly present, 
e.g., as revealed in the BT09 MIR extinction maps, but to gauge the effects of these on the GO emission requires 
higher resolution molecular line maps of the clouds. Eor r ^ 1, Equation ([^ can be simplified to express Tj^ by: 
Tjy = {TB,jyk/[hu])[f{Te^) — /(Tbg)]“^. With this simplification, for an observed voxel Tb,i^ and an assumed Tex, the 
optically thin ^^GO column density per voxel is given by combining Equations Q and ([^: 


dNisco 

dv 


= 1.251 X 10 


14 


Q 


rot 




cm 


-2 


/(Itx) - /(Ibg) 1 - exp{-hu/[kTeA km s-i' 


(3) 


Early studies of IRDCs estimated typical gas kinetic temperatures Tgas ~ 20 K ([Carey et al.|1998 2000 Pillai et al 


2006 ). IRDC E was estimated to have a temperature of 19 K based on NH 3 (1,1) and ( 2 , 2 j VLA observations (Uevme 


2009). However, as discussed below, CO excitation temperatures appear to be significantly lower. 


In our previous study of IRDC H, we used IR AM 30m observations o f C^^O( 2 -l) and (1-0) emission from around 
the IRDC filament to estimate a mean Tex ^ 7 K ( Hernandez et al.|| 2011 [ hereafter H ll). Here, for a uniform analysis 


of the 10 IRDCs, we now use ^^CO Tex estimates from Roman-Dnval et al. ( | 2010 [ hereafter RDIO). In their study 
of 580 molecular clouds, brightness temperatures from CO(1-0) emission line data (Univ. of Massachusetts-Stony 
Brook (UMSB) Galactic Plane Survey), were used to derive proxy ^^CO excitation temperatures, assuming ^^CO 
emission was optically thick and that and ^^CO excitation temperatures are equal. Ultimately, RDIO cited a 

mean excitation temperature for all their molecular clouds based on all clou d voxels above da^ms; The RDIO clouds 
were extracted from the GRS data using a modified version of CLUMPE IND (iWilliams et al.|199 4l), which allowed for 
varying thresholds (i.e., contour increment and minimum brightness, see|Rathborne et al.|(|2009p tor details). 

Eight of our IRDCs overlapped with at least one of their molecular clouds, and in these cases we adopted Tex from 
the RDIO value from the overlapping cloud(s). Eor the remaining two IRDCs (D and E) we set Tex = 7 K, similar to 
the mean values of 7.2 K of Hll and 6.32 K of RDIO. Our adopted values of Tex are listed in TableEor our column 
density estimates that assume optically thin emission, these Tex values are assumed to be constan t throughout 

the cloud. These temperatures are slightly lower (by a few K) than those u sed in previous stud ies (e.g., Simon et al. 


2001 2006b who assumed a fixed value of 10 K). However, the results from Heyer et al. (2009) indicate that GO gas 
throughout GMGs is mostly sub-thermally excited. Note that, in this optically thin limit, varying Tex from 5 K to 
10 K would change the derived column density of the cloud by only ^ 20%. 

Eor the opacity-corrected case, the use of a single mean excitation temperature of relatively low value (^ 7 K) can lead 
to non-physical results in certain regions of the cloud. Equationj^implies = — In {1 — [{TB^yk/[hi'])[f{TQ^) — /(Tbg)]“^]). 
Thus, for a given observed T^^jy, will become undefined if Tex < Tex,crit (he., when [(Tg^jy/c/[hz/])[/(Tex)—/(Tbg)]“^] > 

1). Eor example, a voxel with an observed brightness temperature of Tb,i^ = 5 K will have a numerically undefined 
opacity at Tex,crit < 8-2 K. Hll showed modest temperature variations were present in IRDC H, with a peak temper¬ 
ature of Tex ^ 10 K within the densest clumps and Tex ^ 7 K in the more diffuse gas within the filament envelope 
(see Hll, Eig. 1). Hence, we expect that the voxels containing the largest brightness temperatures have excitation 
temperatures which are a few K larger than the constant Tex values adopted from RDIO. 

To estimate the opacity-corrected column density, we first apply our adopted Tex to estimate r^y in each voxel. Then, 
for each voxel with an undefined r^y, we specified a new excitation temperature of Tex,crit + 1 K, given the voxel’s 
observed T^,zy. This revised excitation temperature allows us to estimate a real solution for dNi^co/dv throughout 
the cloud. After considering a range of possible temperature offsets from Tex,crit, we estimate that the uncertainty in 
dNisco/dv in an individual voxel by this method is at a level of ^ 20 %. 

Ideally, excitation temperatures would be estimated locally from J = 1 —0 observations of the clouds, assuming 

that the line is optically thick. However, there are currently no other surveys that match both the resolution 
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and spatial extent of the GRS. The widely used C olumbia-CfA survey covers the whole Galactic plane, but with 


a low angular resolution of 8' (Dame et ah 2001). The RDIO temperatures measurements are based on the UMSB 


survey, which has a 44" spatial resolution with 3" sampling. 

For the clouds defined by SE extraction and with R = 30 pc, we find that on average 0.7% of the cloud voxels 
require higher temperatures than those listed in Table Cloud D has the highest percentage, 1.5%, and for these 
voxels Tex,crit peaks at 16.4 K with a mean Tex,crit = 8-3 K. For the clouds defined by CE extraction and Ra, we find 
that on average 10% of voxels require higher excitation temperatures. Here, Cloud J has the highest percentage, 27%, 
with a peak Tex,crit of 25.0 K and a mean Tex,crit = 9.2 K. 

The ^^CO-based mass p er voxel, dM was then calculated by assuming a ni 2 co/^i 3 CO = 54 ( Milam et al.||2005 ) and 
^i 2 Co/^H 2 = 2.0 X 10“^ (Lacy et al.||l994). Hence, the assumed abundance of ^^CO to H 2 is 3.70 x 10“^. The mass 
located at a given voxel is then given by 


dM = 1.45 X 10 


-4 dNisco 


dv 


r d y Al Ab 

1013 cm-2 0.212 kms-i J 


Mo, 


(4) 


where d is the cloud distance, Al and Ab are the angular sizes of the GRS pixels, and assuming a mass per H nucleus 
of /iH = 2.34 X 10“^^ g. The total ^^CO-derived cloud mass, M, is simply the total mass of all cloud voxels within 
radius R and velocity range vq ± 15 kms“^. For clouds defined by SE, any pixels with total integrated intensities 
below 5crrms were omitted from further analysis. 

We estimate an uncertainty of 20% in Tex due to its intrinsic variation, which in addition to uncertainties due to 
intrinsic abundance variations, leads to ^ 30% uncertainty in the mass surface density (E). We thus estimate ^50% 
random errors in M, after accounting for these uncertainties in E and the cloud kinematic distance estimates (assumed 
to be ^ 20%). However, we also anticipate that there could be global systematic uncertainties in M (of the whole 
cloud sample) of up to a factor of 2, given the uncertainties in overall absolute ^^CO abundance. 


3.3. Cloud Kinematics 

We used co-added ^^CO column density velocity distributions (e.g.. Fig. Right) to determine the mean velocity, 
vq of each extracted cloud. The velocity dispersion was estimated using two standard methods: 1) the rms ID velocity 
dispersion, a; 2) the width of a fitted Gaussian profile, (Jg- We estimate an mean uncertainty in a of 10%. Additionally, 
we used these Gaussian fitted profiles to estimate a Gaussian profile mass, M^, of each cloud. To visualize how vq and 
(j vary throughout the cloud, we show mass-weighted first and second moment maps of each GMG (e.g.. Fig. |^. 

Using the first moment maps, we derive the velocity gradients in each spatial direction. For example, the longitudinal 
velocity gradient, dvojdl^ was derived by first estimating the mean velocity at each longitudinal position, via a mass- 
weighted sum along the perpendicular (b) direction, then finding the best (mass-weighted) linear fit to these velocities. 
This method was repeated for dvo/db. The magnitude of the total linear velocity gradient across the cloud, dvo/ds^ 
and its position angle direction, were then calculated. Note, we choose to use mass-weighted velocity gradients 
to prevent the results being unduly affected by tenuous wisps of cloud material. Also, in the context of interpreting 
velocity gradients as due to rotation and thus measuring rotational energies of the cloud (below), this mass-weighted 
gradient is the appropriate one to use. 

Many studies of the d ynamics of molecular clouds have interpreted total GMG linear velocity gradients as due to 


solid body rotation (e.g., Philhps|1999 Rosolowsky et al.|2QQ3 Imara fc Blitz|2QlT ). However, the possibility remains 


that the identified single “cloud” actually consists of spatially independent structures. For a cloud undergoing solid 
body rotation, the line of sight velocity gradient in the plane of the sky, dvo/ds, is equal to the projected angular 
velocity, i.e., Uq = dvo/ds. The true angular velocity is U = Uo/sim, where i is the angle between the rotation axis 
and the line of sight. 

The position angle of the projected rotation axis of the cloud contains information that may constrain theories of 
GMG formation and evolution. For example, if a GMG forms rapidly from atomic gas in the Galactic p lane, then 
the GMG rotation is expected to be prograde with respect to Galactic rotation (e.g.. Tasker & Tan 2009). If strong 
gravitational encounters and collisions are frequent between gravitationally bound GlVlCs (Tan 2000), then a more 
random set of orientations of the positional angles of projected rotation axes are expected, including both pro- and 
retrograde rotating clouds. For clouds observed in the Galactic plane, —90° < Oy < +90° represents retrograde rotation 
and —90° < Oy < —180° and 90° < Oy < 180° represent prograde rotation. 

We then also estimated the projected moment of inertia, Jq = ^dMis‘1^ of each cloud using the sky projected 
rotation axis, defined by Oy and the cloud center-of-mass coordinate, where Si is the shortest distance to the rotation 
axis. This allows estimation of the projected rotational energy of the cloud, Trot,o = (l/2)/o^o- 


4. RESULTS 

4.1. Molecular Cloud Physical Properties 

Table [^presents the derived physical properties for all 10 clouds. The global properties of each cloud were evaluated 
for data extracted out to seven different radii, shown on different lines in the table: 5, 10, 20, and 30 pc using SE; 
and out to mass weighted (Rm), areal {Ra)^ and half-mass {R 1 / 2 ) radii using GE. Each cloud definition case is noted 
in column 1. The center-of-mass coordinates of each cloud are listed in columns 2 and 3, evaluated for the optically 
thin and opacity-corrected cases. We note that the GE cloud definitions with Rm, Pa and R 1/2 share the same 
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center. Mean velocities are reported in column 4. Columns 5 and 6 list the separations between the CMC centers 
and the IRDC centers in plane of sky distance and velocity offset, respectively, again for both optically thin and 
opacity-corrected methods. Column 7 lists the cloud radii (with Rm, Ra and R 1/2 cases having two values for the 
optical thin and opacity-corrected cases). Column 8 lists the three cloud mass estimates: optically thin mass (M), 
opacity-corrected mass and the mass estimated from a gaussian fit of the co-added opacity corrected column 

density spectra (Mg)- Column 9 lists the three velocity dispersions measured: the standard deviations of the cloud’s 
optically thin and opacity corrected column density spectra (a, and the dispersion estimated from the gaussian 
fitted opacity corrected column density spectra {(tg)- Column 10 presents the two velocity gradient magnitudes, one 
estimated from the optically thin column density weighted first-order moment map (dvo/ds) and one estimated from 
the opacity-corrected column density weighted first-order moment map {dvo^r/ds). Column 11 presents the position 
angles of the projected cloud rotation axes estimated from each cloud’s angular momentum for both optically thin and 
opacity corrected masses {Oy and The projected moments of inertia estimated using the optically thin mass (Jq) 

and opacity corrected mass (/o,r) ^-^e presented in column 12 . Column 13 shows projected rotational energies, £^rot,o 
and ^rot,o,r 5 and column 14 shows the ratios of these to the gravitational energy, i.e., (3 and Finally, column 15 
shows logio of the virial parameters of the clouds using the three measures of ma ss, i.e., <avir,G- 

A visual overview of each GMC-scale region is given in Figure and Figures |A1||A9[ The left side of the figures 
display a series of images of the molecular gas out to a radius of 30 pc. The first row displays the GLIMPSE 8 /im image 
and the ^^CO mass surface density map . Ei.sco, generated using the velocity interval given solely by the IRDG. The 
IRDC elliptical boundary is also shown ( Simon et al.||2QQ6a ). The subsequent rows display the Eisco maps assuming 
our defined velocity interval of vq ± 15 kms“^, the column-density weighted mean velocity (first moment), the linear 
momentum map derived from each pixel’s column density and velocity with respect to the cloud center of mass, and 
the velocity dispersion map (second moment) for the SE cloud {Left column) and CE cloud {Right column). The right 
side of the figures display the optically thin and opacity-corrected co-added column density velocity spectra at various 
extraction radii for the SE and CE clouds and the SE IRDC. 

Eour of the clouds (D, E, G, and H) were found to contain multiple velocity compone nts, based on the co-added 
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CO column density spectra of their SE 60 pc diameter GMCs (e.g.. Cloud D; Eig. A3). We used a series of 


integrated intensity maps to isolate the velocity range of the molecular gas associated with the IRDC (Eig. |4]). In 


general, these velocity ranges span ^ 10 kms ^ (Cloud D: 84.1 — 93.7 km s Cloud E: 54.3 — 65.5 km s Cloud 
G: 39.7 — 51.8 km s“^; Cloud H: 39.7 — 49.9 km s“^). To investigate how the cloud’s derived physical properties were 
affected by isolating only these velocity ranges, we include these four cases as separate entries in Table (noted with 
“-s”). Eor all subsequent analyses, the “-s” values for these four clouds are used in place of those derived using the 
full velocity range of vq ± 15 km s“^. 

4.2. IRDC and CMC Morphologies 
4.2.1. Description of Individual IRDCs/GMCs 

Cloud A: We estimate a cloud areal radius of Ra — 27 pc and a mass of Mg,a — 2.7 x 10 ^ Mq. The GRS ^^CO GMC 
study by RDIO identified five molecular clouds within the extent of Cloud A, with a total ^^CO mass of 2.1 x 10^ Mq 
(note , we scaled all ^^CO masses of RDIO by a factor 0.49 to convert to our assumed ^^CO/H 2 abundance ratio; see 
^3.2). IRDC A is located about 20 pc north of its parent GMC’s highest E clump, whic h appears bright at 8 fim and 
is thus likely to be star-forming. The WISE H II region catalog ( Anderson et al.||2014 ) lists the W39 H II region to 
be within ^ 3' of this clump’s center, and with a similar LSR velocity (i;lsr = 65.5 km s“^). Another H II region lies 
in a smaller dense clump at I = 19.07°, b = —0.28° with velocity and distance consistent with those of Cloud A. The 
isolated, simple ^^CO emission profiles, along with the spatial connection between the IRDC and the H II regions, 
demonstrated through CE, indicates that GMC A contains clumps at a variety of stages of star formation activity. 

Cloud B: We estimate a radius of Ra — 12 pc and a mass of Mg,a — 3.3 x 10^ Mq. RDIO estimated a ^^CO mass 
of 2.4 X 10^ Mq and a smaller projected areal radius of 6.3 pc. The IRDC is located within the densest region of the 
GMC, which extends ^ 10 pc to the south-east of the IRDC and accounts for most of the cloud’s mass. Spatially, 
Cloud B is well defined and completely isolated at an extraction radius of 20 pc in CE. Although the cloud is isolated 
in velocity space, the column density profile indicates a low level of ^^CO emission at neighboring velocities, both lower 
and higher. Although the GLIMPSE 8/im emission map indicates a large bright emission structure to the west of the 
cloud, this feature is not at velocities consistent with Cloud B (i;LSR ^ 50kms“^). 

Cloud C: We estimate a GMC radius of Ra — 29 pc and a mass of Mg,a — 3.5 x 10^ Mq. RDIO derived a radius 
of 20.9 pc and a mass of ^ 3 x 10 ^ Mq. IRDC C is located within the dense central region of its parent GMC. 
The integrated intensity maps indicate dense and clumpy structure throughout the cloud (Eig. |A 2 [ ). The velocity 
distribution shows a broad centrally- peaked profile , but w ith potential substructure. Additionally, there is a low level 
of ^^CO emission at higher velocities. Butler et al. (2014) studied the IRDC C using near- and mid-infrared extinction 
mapping to probe the high-dynamic range of its mass surface density (down to Ay ^ 1 mag), finding an IRDC mass 
of - 7 X 10^ Mq. This is one of the most massive, hi gh-E IRDCs known, with the potential to form a massive star 
cluster. Using the same IRDC elliptical boundary from Simon et al. (2006a), we estimate a ^^CO mass of 5.7 x 10 ^ Mq. 
Our slightly lower mass estimate may be the result of abundance uncertainties and/or CO depletion. 

Cloud D: Eor the GMC at 'Ulsr = 84.1 — 93.7 km s“^, we estimate a radius of Ra ^ 21 pc and a mass of Mg,a — 
9.3 X 10^ Mq. RDIO identified two clouds within the Cloud D-s region with a total ^^CO mass of ^ 5.2 x 10^ Mq. 
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Fig. 3.— Cloud A: Left panel maps: Top row: Left: GLIMPSE 8/i,m image. Right: mass surface density (S) map over the 

velocity range of the IRDC emission profile v\sr = 54.6 — 77.1 kms“^ (ellipse denotes IRDC boundary). 2nd row: Left: SE ^^CO-derived 
S map (circles denote 5, 10, 20 pc boundaries). Right: CE ^^CO-derived S map. 3rd row: Left: SE mean velocity (first moment) 

map. Right: CE ^^CO mean velocity (first moment) map relative to vq defined from the IRDC. The arrows show the direction of total 
(mass-weighted) line-of-sight velocity gradients, dvo/ds. An arrow is shown for each extraction radius (defining the length of the arrow) 
and is centered on either the extracted cloud’s center of mass (for SE cases) or the R = 30 pc cloud’s center of mass (for CE cases). 4th 
row: Left: SE ^^CO-derived (optically thin) linear momentum map. Right: CE ^^CO-derived (optically thin) linear momentum map. The 
arrows show the position angle, Oy, of the cloud’s rotation axis, i.e., orthogonal to the direction of the velocity gradient, for each extracted 
cloud (with arrow length showing cloud radius) 5th row: Left: SE ^^CO velocity dispersion (second moment; optically thin) map. Right: 
CE ^^CO velocity dispersion (second moment; optically thin) map. All maps in Rows 2-5 cover the velocity interval of no ± 15 kms“^, 
with no defined by the emission of the IRDC. Omitted pixels, those that fall below the 5crrms level, are shown in gray. Right panel veloeity 
distributions: In each panel, the column density distributions with velocity are derived assuming optically thin conditions (blue solid line) 
or with opacity correction (red solid line) from the CO(1-0) co-added spectra. Gaussian fits are shown by the dotted blue and red lines. 
In all panels, a and ac are noted in the upper left corner of the plot. Left eolumn: Top to bottom: SE out to R = 5,10, 20, 30 pc. Right 
eolumn: Top to bottom: SE of IRDC ellipse, CE out to Rm? Rif2- 
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Fig. 4.— channel maps of the multi-component clouds, top left: Cloud D; top right: Cloud F; bottom left Cloud G; bottom right: 

Cloud H. The velocity interval of each panel is shown in the upper-right corner of each map. 


The IRDC is located in the dense central region of the GMC. However, the spectra show four cloud components 

within vq ± 15 kms“^. Although the emission corresponding to the IRDC is evident within a radius of 5 pc, 

a lower-velocity gas component begins to dominate on scales out to about 20 pc. The position-velocity map (Fig. 

shows there is a significant amount of emission at lower velocities, t’LSR = 63.5-84.1 kms“^, that extends 

largely over the northern region of the cloud. In the same velocity interval of the IRDC, '^^lsr = 84.1 — 93.7 km s“^, a 
substantial molecular gas clump contains the largest integrated peak intensity 50 K kms“^) at the western edge 
of the 30 pc extraction radius. Figure [^suggests that this emission is only tenuously connected to the ^^CO emission 
associated with the main CMC. 

Cloud E: We estimate a GMC radius of Ra — 28 pc and a mass of Mq^a — 3.1 x 10^ Mq. No clouds within the 
Cloud E spatial and velocity ranges were listed in the RDIO catalog. The integrated emission map indicates that the 
GMC is highly substructured. The IRDC is contained within a small dense clump central to its parent GMC, but 
more massive clumps are seen on scales beyond 5 pc. CE finds an extended region of clumpy ^^CO gas in the cloud’s 
south-west region, beginning at 10 pc and extending beyond the 30 pc boundary. 

Cloud F: Eor the cloud in the velocity range t’LSR = 64.3 — 65.6 kms~^, we estimate a r adius of Ra ^ 21 pc and 
a mas s of Mq^a — 7.8 x 10^ Mq, assuming a distance of 3.7 kpc (Rathborne et al. 2006). Note, Knrayama et al. 


( 2011 ) estimated a parallax distance 1.56 ±0.12 kpc, however this result has been called into question by Foster et al. 


(2012). RDIO found two GMCs associated with Cloud E: one cloud matches the location and velocity interval of Cloud 
E’s western massive clump, having a mass of 6.9 x 10^ Mq and radius of 12.3 pc. The other molecular cloud is 
associated with the smaller eastern clump, for which they estimate a mass of ~ 2.2 x 10^ Mq. Our analysis recovers 
this larger total mass, ^ 10^ Mq, if we include all the ^^CO emission within vq ± 15 km s“^. The larger molecula r cloud 
complex sur rounding Cloud E was recently cataloged within the Giant Molecular Eilament (GME) study of Ragan 


et al. (2014). They found that their filament GME38.1-32.4 is actually the projection of two overlapping filamentary 
structures. Cloud E is within the further filament (3.0-3.7 kpc), which has a total length of 232 pc. The GME contains 
multiple infrared dark regions (identified in both GLIMPSE Sfim and HiCAL 250/im images) and several dense gas 
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clumps (identified by the Wienen et al. ( 2012 ) NH 3 and Bolocam Galactic Plane surveys), demonstrating that overall 
it is in an early phase of star formation. However, in Cloud F alone, there are two H II regions within the two densest 
parts of the cloud wit h consistent LSR velocities and distances, at I = 34.256,6 = 0.136 and I = 34.821,6 = 0.351 
(Anderson et al.|2014|). While IRDC F, which is also highly filamentary, contains some of the highest column densities 
in t he CMC, the highest S clump is actually another feature extending west of the IRDC between 5 and 25 pc (Fig. 
A5). The velocity distribution of the GMC for material out to 30 pc indicates that there are four cloud components 
ithin our defined velocity interval uo±15 km s“^. Figurej^shows the PV maps of the cloud out to 30 pc and indicates 
that the majority of the ^^CO emission resides within the velocity interval of the IRDC (ulsr = 54.3 — 65.6 kms“^) 
and is connected to the weaker molecular emission at surrounding velocities. Additionally, the massive clump and 
IRDC are kinematically associated with another dense clump on the eastern edge of the cloud at ^ 25 pc. The cloud 
defined by CE out to a radius of 30 pc indicates that the regions above and below the IRDCs parent cloud are filled 
with clumpy molecular material related to the other three emission profiles. 

Cloud G: For the material in the velocity range of v\sr = 39.7 — 51.8 km s~^ we estimate a cloud radius of Ra — 22 pc 
and a mass of Mq^a — 8-0 x 10^ ATq. RDIO identified two clouds within the Cloud G boundary with a total ^^CO mass 
of ^ 5.7 X 10^ Mq. At small cloud radii, < 10 pc. Cloud G’s velocity distribution is reasonably well described by a 
single peak. However, at larger radii we clearly see three separate components within uq ± 15 km s“^. The IRDC, and 
its parent GMC, are confined to the central velocity interval oiv = 39.7 — 51.8 km s“^. In the GMC defined by SE, the 
IRDC lies at the northern edge of the cloud’s densest region, which extends south ^ 18 pc. Eor the cloud defined by 
CE, this same cloud region breaks up into three smaller clumps. The lower velocity component of the cloud contains 
a dense m olecu lar clump at I ^ 35.2, 6 —0.7, a region which also contains a high level of mid-infrared emission at 

8 /im (Eig. |A 6 [ ). The WISE H II region database lists a candidate H II region near the center of the low-velocity clump 
with an estimated parallax distance of 2.2 kpc, versus the Cloud G distance of 2.9 kpc. 

Cloud H: Eor gas in the velocity range of uisr = 39.7 — 47.1 km s“^, we estimate a cloud radius of Ra — 19 pc and a 
mass of Mg^a — 5.6 x 10^ Mq. RDIO identified three molecular clouds at the location of Cloud H, with a total mass of 
^ 2.8 X 10^ Mq. The PV diagrams indicate that there are three separate cloud components within vq ± 15 km s“^ (Eig. 
1^. The IRDC and its parent GMC are confined to the velocity interval of u = 39.7 — 47.1 km s“^, in which the IRDC 
IS the main structure of a larger diffuse cloud extending to the west. There are additional dense gas clumps at the 
neighboring velocities: a large coherent clump about 25 pc to the south at u = 32.5 — 39.7 km s~^ (part of Cloud G), 
and a high concentration of dense, clumpy gas along the entire northern region at velo citie s of v = 47.1 — 62.6 km s“^. 
These higher velocity clumps overlap with a region of bright 8/im dust emission ( Eig. |A7|). The WISE H II database 


lists five H II regions overlapping these dense clumps with similar LSR velocities (Anderson et al. 2Q14[ ). 


Cloud I: We find that GMC I has a radius of Ra — 12 pc and a mass of Mq^a — 2.6 x 10"^ Mq. However, RDIO 
identified two clouds that match the location of Cloud I, and estimated the much larger mass of ^ 2.4 x 10^ Mq. 
This difference is due to the surro unding diffuse ^^CO eui ission, which is included by their cloud finding method that 
uses the CLUMPEIND algorithm (jRathborne et al.||2009|). IRDC I lies within the densest region of the GMC, which 
extends roughly 7 pc north of the IRDC. 'I'he cloud's SE column density spectra indicates a symmetric and well 
isolated emission profile out to radii of ^ 10 pc. However, at larger cloud radii, diffuse ^^CO gas at higher velocities is 
seen. However, the spectra of the cloud defined by CE indicates that this high-velocity emission is at intensities less 
than the 5 cr threshold, suggesting that the GMCs is surrounded by a reservoir of low-density molecular gas. Eor the 
cloud defined by CE, the central dense region is broken down into three dense clumps and the GLIMPSE 8 /im image 
indicates that there are multiple bright dust regions along these dense ^^CO clumps. However, only one H II region 
is identified within the WISE H II catalog, with a slightly lower LSR velocity of ^ 38 km s“^. 

Cloud J: We find that Cloud J has a radius of Ra — 13 pc and a mass of Mg,a — 2.0 x 10^ Mq. However, Cloud J 


was included in the GME catalog by Ragan et al. ( |2Q14 ) in which they estimate a total cloud length of 79 pc. They 
estimated a total cloud mass of 6.8 x 10"^ Mq (part of this difference results from the larger considered area; part 
from assuming a smaller ni2Co/^H2 ratio of 1.1 x 10“^). Our analysis also finds that at the largest extraction radii. 
Cloud J is a filamentary cloud, extending ^ 60 pc in the CE case. Additionally, the co-added distribution of column 
density with velocity at all cloud radii indicate that C loud J is quite symmetric and well isolated in velocity space. 
IRDC J, i.e., G053.11+00.05 from Simon et al. (2006a), defined the search region, but Cloud J also includes the IRDC 
G053.31+00.00 within its 30 pc rad ius. The 8 /im images of the Cloud J reveal a significant amount of bright dust 

(|2014) associated with a foreground high-mass star-forming region at ulsr ^ 4 kms“^. 
database found on e radio quiet H II region at I = 53.09,6 = 0.12 with an LSR 


emission, which Ragan et al. 

Our search through the WL^E H H 


velocity consistent with Cloud J ( Anderson et al.||2014 ). 


4.2.2. Spatial and Velocity Offsets between IRDCs and GMCs 

The spatial and velocity offsets between the IRDC and GMC center of masses are listed in Table quoted to the 
nearest resolution element of the GRS position-velocity grid. In general, spatial offsets grow as one considers larger 
GMC scales. At the 30 pc aperture scale applied to conn ected extraction, this offset is typically ^2 — 10 pc. All but 
clouds C and I have this GMC center located outside the Simon et al. (2006a) IRDC boundary. 

Velocity offsets are typically relatively small, ^ 1 km/s. This suggests that IRDCs have been formed by contrac¬ 
tion/compression of structures that were already pre-existing within the GMC. Clouds A, B, C, E, I, J have relatively 
simple CE CO morphologies in position velocity space. The remaining clouds exhibit more complex, multi-component 
structures. 




























14 



A IRDC • GMC 

▲ A 
A B 

C 

A D -s 

A E 

A F-s 
A G -s 
A H-s 

A I 
A J 


Fig. 5.— Velocity gradient dvo/ds versus radius R of the clouds over which the gradient is evaluated. Each panel represents a different 
cloud definition: (a) Top left: simple extraction, optically thin column density derivation; (b) Bottom left: simple extraction, opacity- 
corrected column densities; (c) Top right: connected extraction, optically thin column density derivation; (d) Bottom right: connected 
extraction, opacity-corrected column densities, n each panel, the 10 IRDCs are shown by the colored filled triangles as defined in the upper 
right-hand corner. For Clouds D, F, G, and H, these estimates are based on th e single component emission profile. GMCs are shown 
with filled circles. The best-fit power-law relation from [McKee fc: Ostriker ( |2007| {dvo/ds ~ (7pc(s/pc)'^'^“^, with values of aa — 0.5 and 
cTpc — 0.72 kms“^) is shown by the dotted lines. The global (IRDCs and GMCs, excluding IRDC J) best-fit power-law relation for each 
panel is shown by the solid lines. 


4.3.1. 


4.3. Kinematics and Dynamics 
Velocity Gradients and Implications for Turbulence and Rotation 


The GMCs defined by CE have a mean velocity gradient of dvo^r/ds ^ 0.14 km s ^ pc A similar value of 


0.13 kms ^ pc ^ was found by |lmara fc Blitz (2011) for five GMCs. The 

velocity gradients that ar e smaller than these values. Jackson et al. (2010) find dv/ds < 0.09 km s " p c 


100-pc-scale filaments have global 
-1 ^.-1 in “Nessie”, 


Battersby & Bally (|2012 ) find dv/ds < 0.05 kms ^ pc ^ in their 80-pc long cloud, and Ragan et al. (2014) measure 

r\rty / O on Ervi o -tA Z' • -L o o on OTmT*orm /of "Flm V -Rlom/on-fo in -flimT* oomnl/o 


di;/d5 ^ 0.06 km s ^ pc ^ as an average of the 7 filaments in their sample. 

For the elliptical IRDC regions also defined by CE, we find a mean velocity gradie nt of dvo.T/ds 0.26 k m s“^ 


larger than that of the GMCs. For comparison, on smaller scales within IRDC H, Henshaw et al. (2014) find global 
velocity gradients of 0.08, 0.07, 0.30 km s“^ pc“^ in several sub-filaments measured over 


pc 


-I 


on centroid velocities of the dense gas tracer N2H +(1 — 0). Larger local gradients 
present on sub-parsec scales. Ragan et al. ([2012) found velocity gradients of 2.4 and 2.1 kms“^ 

regions of IRDCs G5.85-0.23 and C24. 05-0.22, based on observation s of NH^ (1, !)• __ 

Observed velocity gradients may be a signature of turbulence (e.g., Burkert & Bodenheimer|2000 McKee & Ostriker 

TT ,1. .-1,1 /T -r-T’rr-i \ i , • -i n )) A K JKj I ^ I ^ ' n W I 


^2 parsec length scales, based 
1.5 — 2.5 kms“^ pc“^ can be 
within sub-pc 


^pc 


2007). For a turbulent line-width-size ( LWS) relation of the form a = crpc(s/l pc)"^^, where values of aa- — 0.5 and 
CTpc — 0.72 kms“^ have been reported (McKee & Ostriker 2007), we expect velocity gradients measured over some 
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Fig. 6. — Velocity dispersion a versus radius R of the clouds over which the dispersion is evaluated. Each panel represents a different 
cloud definition: (a) Top left: simple extraction, optically thin column density derivation; (b) Middle left: simple extraction, opacity- 
corrected column densities; (e) Bottom left: simple extraction, opacity-corrected column densities fitted by a Gaussian profile; (d) Top 
right: connected extraction, optically thin column density derivation; (e) Middle right: connected extraction, opacity-corrected column 
densities, (f) Bottom right: connected extraction, opacity-corrected column densities fitted by a Gaussian profile. In each panel, the 10 
IRDGs are shown by the colored filled triangles as defined in the upper right-hand corner. For Glouds D, F, G, and H, these estimates 
are based on th e single component emission profile. GMGs are shown with filled circles. The best-fit power-law relation from |McKee &:| 
|Ostriker| (|2007| (cr = (7pc(s/l , with values of aa — 0.5 and fjpc — 0.72 km s“^) is shown by the dotted lines. The global (iHDGs and 

VilVUJsj best-til power-law relation for each panel is shown by the solid lines. In panel (f), the best-fit power-law relations are also shown 
for the IRDGs (dot-dashed line) and the GMGs (dashed line). 
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r-l 


scale s to vary as dvo/ds crpc(5/pc)^ 

Figure [^displays our estimates of dvo/ds as a function of cloud size for both IRDCs and GMCs (we note that the 
small size of IRDC J makes it difficult for velocity gradients to be properly resolved with GRS data, so we exclude 
it from the following analysis). Fitting the data for clouds defined by CE and including opacity-corrected column 
densities (lower-right panel) we find ao- = 0.39 ± 0.14 and a^c = 0.54 ± 0.19kms“^. If each cloud is considered 
individually, has mean and median values of 0.45 and 0.51, respectively, and a standard deviation of 0.59. By 
inspection of the other panels in Fig. we see that the results for the slope and normalization of the power law relation 
are not too much affected by the choice of cloud definition and mass measurement, but the GE results show smaller 
dispersions. 

We also examine directly the observed velocity dispersions as a function of size scale (effective radius) of the clouds. 
Eigure displays estimates of cr as a function of cloud size for both IRDGs and GMCs. Now we see a larger effect 
on the velocity dispersion versus size relation depending on how the cloud is identified (i.e., SE or CE), how its mass 
is measured (i.e., whether or not to carry out an opacity correction for the ^^CO-derived mass), and how to derive 
velocity dispersion (i.e., via a direct estimate or via a fitted Gaussian). The CE, r-corrected, Gaussian-fitted results 
yield a relatively tight power-law, with the smallest dispersion. 

Eitting the data for all clouds defined by CE and including opacity-corrected column densities (lower-right panel) 
we find = 0.25 ± 0.06 and a^c = 1-07 ± 0.15 km s“^. If each cloud is considered individually, has a mean and 
median value of 0.17 and 0.16, respectively, and a standard deviation of 0.10. If we fit IRDCs and GMCs separately. 


we find an- = 0.27 ±0.11 and a. 


pc 


1.17 ± 0.16 km s ^ for IRDCs and a^r = 0.58 ± 0.13 and a, 


pc 


0.43 ±0.15 km s- 


for GMCs. Eor virialized, constant mass su rface density clouds, where a oc (GM/R)^/^ oc y;i/2/^i/2 (|Larson| 


1981 


McKee fc Ostrih^|2007 McKee et al. 12010 ) one expects a^j = 0.5 and this is consistent with the GMC-only results 


However, lor IRDCs alone or lor the full sample including IRDCs, we see a shallower relation, which could still be 
consistent with expectations from virial equilibrium si nce IR DCs have higher mass surface densities. We will discuss 
the virialization of the clouds in more detail below in §4.3.3 


Observed velocity gradients w ithin molecular clouds can also be interpreted as due to cloud rotation (e.g., Goodman 
et al. II 1993 Imara fc Blitz||2Qll ). With such an interpretation, we assess the projected rotational energy of the clouds 
Erot,o = (l/2)/or^o5 discussed in §3.3. The derived values of Iq an d Erot.n are listed in Tabl e 2. 

We next compare rotational to gravitational energies, 
gravitational energy of an ellipsoidal cloud is given by 


Eollowing (Bertoldi & McKee 1992 , hereafter BM92), the 


w 3 GM2 

W = -aia2—^, 
5 R 


(5) 


where parameter ai describes the effects of the internal density distribution of the cloud, where for a power-law 
distribution p oc ai = (1 — /cp/3)/(l — 2kp/b). Here, for both GMCs and IRDCs, we adopt kp = 1, such that 

ai = 10/9, based on the IRDC density profile study by BT12. Somewhat steeper density profiles with kp 1.5 are 
derived by BT12 if envelope subtraction is carried out, but this would only change ai by ^ 10%. The parameter a 2 
accounts for the effect for the cloud’s ellipticity. Eollowing BM92 and our ellipsoidal cloud virial analysis from HTll, 
an ellipsoidal cloud with radius R normal to the axis of symmetry and size 2Z along the axis has an aspect ratio of 
y = Z/R^ and semi-major and semi-minor axes of rmax and rmin, respectively. Since we do not know the inclination 
angles of the IRDCs or GMCs, we assume that R = rmin, Z = rmax, and that the observed cloud radius is equivalent to 
the geometric mean radius, Robs = (RminRmax)^^^- These assumptions differ from those of HTll, where we adopted 
fiducial inclination angle values of 60° for the filamentary IRDCs E and H. Now, with a larger sample of both clumpy 
and filamentary IRDCs, we make no assumption regarding their inclination angles. However, as noted in HTll, an 
uncertainty of 15° in inclination would lead to a ^ 15% uncertainty in the value of y. We follow the BM92 definition 
of a 2 , where 

Rrr, arcsinh(2/^ - 


a2 = 


R 


(y2 _ 1)1/2 


where, Rm is the mean value of Robs averaged over all viewing angles. Eor clouds with low aspect ratios {y < 10), 
BM92 give an approximation for the mean observed radius of RmiR ^ i 0.27|1 — ^-2|0.8i^o.3i^ accurate to 

within 2%. Since all of the IRDCs in this study have y values below 10, with IRDC B having the largest aspect ratio 
of ^ = 4.67, the BM92 approximation of a 2 is valid for use in our gravitational ener^ estimates. 

We thus evaluate VF, and then the ratio [3 = Rrot/|kF|, which we also list in Tables]^ and We find that the GMCs 
defined by CE and with opacity-corrected column density estimates have mean ^ 0.03^ Eor the elliptical IRDCs 
also defined by CE , we find mean — 0.051. These values are larger tha n the /3 ^ 3 x 10~^ va lues estimated by 
Ragan et al. (2012) for two IRDCs. On much smaller scales, the study by Goodman et al. (1993) found that most 


dense cores within dark clouds have values of (3 < 0.18. The increasing values of [3 that we see on going from CMC to 
IRDC scales may be indicative of slightly more disturbed kinematics within IRDCs. However, it could also be caused 
by systematic changes in mass estimation from ^^CO, e.g., due to increased CO freeze-out within IRDCs. 

4.3.2. “Rotation” Directions 

Eigure[^ displays the distribution of position angles of cloud angular momentum vectors, with respect to the 
direction of Galactic rotation for the cloud definition cases assuming opacity corrected column densities. Our results 
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Fig. 7. — Distribution of relative position angle of the angular momentum vectors of the 10 IRDCs/GMCs assuming opacity corrected 
column densities with respect to the direction of Galactic rotation, The effect of diffe rent cloud defin i tions is illustrated in the different 

panels. In each panel, the black line shows the results from the numerical simulation of|Tasker & Tan| (|2009|) for clouds at t = 250 Myr 
and the appropriate inclination of viewing angle of 0° with respect to the galactic plane, i.e., an m-plane view. 


show that while a cloud’s position angle can vary depending on cloud definition, the overall, global distribution of 
cloud position angles appears consistent with near equal fractions of clouds with prograde and retrograde rotation, 
although based on a relatively small sample of only 10 analyzed clouds. 

Young GMCs, recently formed from the diffuse ISM should possess position angles which are prograde with Galactic 
rotation {Oy < —90°; 6y > 90°). For all of the clouds, measured at various radii for the SE and CE definitions, we find 
that there is a large fraction of clouds with retrograde rotation (—90° < Oy < +90°). We find that ^ 50% of the GE 
clouds analyzed with opacity correct ed column densities have retrograde motion. 

In the Milky Way, Phillips (1999) compiled a database of rotational measurements for 156 individual clouds from 
previously published studies, finding that there is a random distribution between the direction of a cloud’s angular 
velocity and the Galactic plane. However, isolated clouds had orientations towards the Galactic poles. They concluded 
that the rotation of large-scale structures was due to Galactic shear, while the randomized orientation of cloud position 
angles was due to dyna mical and/or magne tic interactions. By estimating the velocity gradients for more than 500 


clouds within the GRS, Koda et al. 


nearly equal. 
in retrograde rotation w i th respect to the g alaxy. 


r |2006D found that the p opulations of prograde and retrograde rotation were 
Rosolowsky et al.| (|2003|) found that 40% of the GMG population in M33 are 
in a subsequent study of the H I envelopes of 45 GMGs from 


Rosolowsky et al. (2003), Imara et al. (2011) found that 53% of these GMCs are in retrograde rotation with respect 


to the local H I. Eurtherniore, they suggest that if the linear gradients of the H I envelopes are due to rotation, then 
62% are in retrograde rotation with respect to M33. 


In comparison to theoretical models, the disk galaxy simulations of Tasker & Tan (2009) found near equal amounts 
of apparent prograde and retrograde clouds 50%) for simulated clouds at t = 250 Myr and a viewing angle of 0° 
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inclination to the galactic plane, i.e., an in-plane view. A substantial fractio n of retrograde clouds re sults if GMCs 
live long enough to undergo mergers and/or collisions with neighboring clouds. |Dobbs & Pringle| (|2013|) performed 3D 
smooth particle hydrodynamic (SPH) simulations of GMCs within a galactic disk, including self-gravit y and stellar 
feedb ack. They found a retrograde GMC population of ^ 40%, 

2011). However, the most massive GMCs tend to be prograde. 


similar to their previous studies (e.g., Dobbs et al 


"Studying the distribution of angular mo mentum v ectors can test scenarios of star formation that involve frequent 
galactic shear driven cloud-cloud collisions (TanpOOO). At the moment, with only 10 clouds in our pilot sample, we are 
not yet in a position to make detailed statistical comparison with numerical simulations. We defer such an analysis, 
using larger samples of Galactic GMCs, to a future paper. 


4.3.3. Cloud Virialization 

Various mechanisms, such as turbulence, magnetic fi elds, and thermal gas pres sure, may be responsible for supporting 


molecular clouds against gravitational collapse (e.g., McKee fc Ostriker 2007). The dimensionless virial parameter 


Q^vir, describes a cloud’s dynamical state via ( Bertoldi V McKee||1992p : 


_ 5a'^R Ek 
GM ~ \Eg\ 


(7) 


where a is the ID mass-averaged velocity dispersion, a = \Eg\/{SGM‘^/[ bR]) is the ratio of gravitational energy. Eg 
(assu ming negligible external tides), to that of a uniform sphere, and Ek is the kinetic energy. As discussed above in 
for spherical clouds with power-law density distribution, p oc r~^p 
4, 10/9, 5/4, 5/3 when kp rises as 0, 0.5, 1, 1.5, 2, respectively. 


p.3 


= (1 — kp/3)/{l — 2kplb), i.e, rising as 1, 


Effects of nonspherical geometry for ellipsoidal clumps have also been considered by BM92, and found to be relatively 
small (< 10%) for aspect ratios Z/i? < 6, where Z is the clump radius along the axis of symmetry and R the radius 
normal to this axis. Therefore, here we will ignore these effects of clump elongation. 

Ignoring surface pressure and magnetic terms, a cloud in virial equilibrium has \Eg\ = 2Ek^ so that o^vir = cl- Such 
a cloud in free-fall has 2a from below as time progresses. In principle one can imagine a cloud starting with 

very small levels of internal motion, so that 1, and then crossing through the condition = a as part of its 

free-fall collapse. However, such a tranquil and synchronized initial condition would seem to be difficult to achieve in 
reality if GMGs are accumulating from a turbulent ISM, including potential agglomerating via GMG collisions. 

Although it is difficult to differentiate infall motions from vir ialized turbulent m otions, in cases where infall rates 
have been claimed towards dense clumps, including IRDGs, (e.g., Peretto et al.|2Q13 ), the infall velocities are relatively 
small compar ed to virial motion s (equivalent to the mass doubling times being lo ng (^ 10 x) compared to the free- 


fall time; see Tan et al.1 (2014), for a review). On the larger scales of GMGs, Zuckerman & Evans (1974) and 


Krumholz V 'ran~( 2007p J rave argued against fast rates of global collapse and infall, since overall star formation rates 
are much smaller. However, converging flows associat ed with the collision of independent, self-gravitating clouds may 


V J.J.WVVO CAjO OW v^xcxuv 

be relatively common ( Tan||2QQQ Tasker fc Tan||2QQ9 ), and the kinematic signatures of such flows would be difficult to 
distinguish from global infalT 

The effect of surface pressure on the cloud increases the velocity dispersion of the virial equilibrium state and thus 
raises the value of its virial parameter. Eor example, for a cloud with kp = 1.5 and negligible surface pressure we 
have ayir = cl = 5/4. However, if the cloud is embedded in an ambient medium wit h the cloud’s surface pressure 
equal to the ambient pressure, as in the Turbulent Gore model of McKee & Tan (2003 hereafter MT03), then a^ij. = 
(5/2)(3 — kp)/[{b — 2kp){kp — 1)] = 15/4. Note as kp ^ 1 from above, a^ir co (a self-similar singular polytropic 
sphere with kp = 1 would have a spatially constant pressure, i.e., zero pressure gradient, so there is no self-consistent 
physical solution; MT03). These results show that virial equilibrium pressure confined singular polytropic spherical 
clouds with kp < 1.5 have virial parameters ^ 1. 

Support by large-scale magnetic fields acts to reduce the velocity dispersion of the virial equilibrium state. Eor 
example, in the fiducial MT03 Turbulent Gore model, the mean velocity dispersion of the virial equilibrium state is 

reduced by a factor ^ 2.8“^/^ = 0.598 if allowing magnetic fields strong enough that the mean Alfven Mach 

number is unity. This implies that the value of the virial parameter of the virial equilibrium state would be reduced 
by a factor of 0.357. 


Magnetic field observations in molecular clouds have been reviewed by Grutcher (2012) and Li et al. (2014). Eield 


strengths increase with density approximately as H oc and are close to values expected for equipartition magnetic 
and turbulent energy densities and also close (within a factor of ^ 2) of the critical field strength that would prevent 
collapse. It thus seems likely that magnetic fields can have significant dynamical effects. Recently, Pillai et al. (2015) 
have found evidence for strong, dynamically important magnetic fields in two IRDGs, based on relatively ordered 
directions of polarized sub-millimeter emission, thought to be caused by dust grains that have aligned with the B- 
fields. Unfortunately for the IRDGs/GMGs in our sample, we do not have detailed information on magnetic field 
strengths. 

Eor the GMGs and IRDGs analyzed in this paper, for each cloud, three virial parameters were estimated. Eirst, 
we estimated assuming optically thin conditions to derive a and M. Second, using the opacity corrected column 
densities, (avir,r is estimated using and M^. Einally, the virial parameter was estimated using the Gaussian-fitted, 
opacity corrected column density spectra, i.e., (avir,G is estimated using gg and Mg- In the cases of clouds D, E, G, 
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Fig. 8. — Comparison of cloud virial parameter, Oivir, with cloud mass. Each panel represents a different cloud definition: (a) Top 
left: simple extraction, optically thin column density derivation; (b) Middle left: simple extraction, opacity-corrected column densities; 
(e) Bottom left: simple extraction, opacity-corrected column densities fitted by a Gaussian profile; (d) Top right: connected extraction, 
optically thin column density derivation; (e) Middle right: connected extraction, opacity-corrected column densities, (f) Bottom right: 
connected extraction, opacity-corrected column densities fitted by a Gaussian profile. In each panel, the 10 IRDGs are shown by the colored 
filled triangles (see legend). For Glouds D, F, G, and H, these estimates are based on the isolated single fitted emission profile. The error 
bar in the top-left of each panel represents the mean uncertainties of ~ 40% in avir and ~ 50% in M. 


and H, the “-s” mass and velocity dispersions were used to estimate the single component virial parameters of the 
clouds. If cloud properties are assessed using the full vq ± 15kms“^ velocity range (including multiple components), 
then larger values of M and especially a are inferred, leading to estimates of that are larger by factors of about 3. 
All molecular cloud virial parameters estimates are listed in Tables and We note that for CE-defined GMCs, the 
virial analysis is carried out using cloud centers that are the center-of-mass of all the CE gas within the 30 pc radius 
region. The cloud masses are those that are estimated from the material that is within the radii Rm^ Ra and 
Eor SE-defined GMCs, the virial analysis uses the extracted cloud radius as the characteristic size of the cloud. 

Eor the same cloud definition cases, the virial parameter was also estimated for each IRDC. Since these clouds are 
elliptical, we used the geometric mean observed radius, R = (Rmaj^min), as an estimate of the cloud’s equivalent 
spherical shape (but note that, as discussed above, effects of asphericity are minor). 

Eigure displays the virial parameters as a function of mass for each cloud definition case. The overall range of 
Q^vir is nearly the same for both IRDGs and GMCs. We estimate a typical random uncertainty in of ^ 40% after 
accounting for uncertainties in the mass surface densities, velocity dispersions, and the cloud kinematic distances. 
However, systematic uncertainties in mass surface density of a factor of ^ 2 would lead to similar levels of uncertainty 
in (avir- Considering the different cloud definition methods, the virial parameters estimated using the CE are lower 
than those estimated using SE. Eor both cloud extraction methods, o^vir is generally smaller for the clouds using 
opacity corrected column densities compared to those assuming optically thin gas conditions. This trend results from 
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Fig. 9.— Dependence of with mass surface density, S. Results for the IRDCs (Left) and GMCs (Right) are shown separately. 

Each panel represents a different cloud definition: (a) Top left: simple extraction, optically thin column density derivation; (b) Bottom left: 
simple extraction, opacity-corrected column densities; (c) Middle right: connected extraction, optically thin column density derivation; 
(d) Middle right: connected extraction, opacity-corrected column densities, (c) Bottom right: connected extraction, optically thin column 
density, Gaussian fitted profile derivation; (d) Bottom right: connected extraction, opacity-corrected column density, Gaussian fitted profile 
derivation. In each panel, the IRDGs/GMGs are shown by the colored filled triangles/circles as defined in the upper right-hand corners. 
For Glouds D, F, G, and H, these estimates are based on the single component emission profiles. The error bar in the upper-left corner 
represents the mean uncertainties of ~ 14% in cr/R^/^ and ~ 3 0% in S. The re sults of RDIO, a total of 329 clouds (M > 10"^ Mq) scaled 
by a factor of 0.49 to reflect our adopted ^^GO abundance (see |Tan et al.|20I3| are shown by the g ray crosses. The dotted line represents 
virialized conditions with avir = 1? he., cr/ijl/2 = (7rG/5)l/2£W (see text). The mean results of 
(black open square). Best-fit power-law relations (cr/R^/^ oc S’^) are shown by the solid lines. 


Solomon et al. 


( 1987t are also shown 


larger masses and slightly narrower line profiles in the opacity-corrected cases. Similarly, the smallest ttvir values were 
estimated in the cases of using the Gaussian htted spectra to measure both the velocity dispersion and mass, cfg and 
Mg- 

Thus, the virial parameter is sensitive to the method used to extract the emission associated with a GMC. 

Overall, the lowest values of o^vir were estimated for the clouds dehned by GE, opacity-corrected column densities, and 
Gaussian htted masses and velocity dispersions. The mean and median values of these estimates of avir,G are 1.12 and 
0.99, respectively, with a standard deviation of 0.51 (similar to the expected size of the random errors). Similarly, for 
the IRDGs under the same dehnitions, the mean and median values of (avir,G are 1.88 and 1.93, respectively, with a 
standard deviation of 1.11. 

We see that IRDGs have virial parameters that are about a factor of two larger than GMGs. Again, this can be 
interpreted as evidence for more disturbed kinematics. Alternatively, it could indicate that ^^GO-derived masses for 
IRDCs have been underestimated, potentially due to CO depletion. 

A population of virialized clouds with similar degrees of pressure conhnement may be expected to have a relatively 
small dispersion in their values of Gvir, although some variation is likely due to variation in internal density structure 
and level of support by large scale magnetic helds. If we denote the mean virial equilibrium value of the virial 
parameters of such a cloud population as Q;vir,eq 5 then from eq. Q we can derive 


a 



1/2 


a 


1/2 

vir,eq 


^ 1/2 


( 8 ) 


(see also Heyer et al. 2009 
GMC samples, we can see i; 


McKee et al. 2010). By examining the observables and E of our IRDC and 

they follow this correlation that is expected if the clouds are virialized. We estimate an 
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Fig. 10.— Dependence of with mass surface density, S, for the all clouds defined by CE, opacity corrected column densities, 

and gaussian fitted velocity dispersions. As in Figure]^ IRDCs are shown by the filled triangles and the GMCs by the filled circles (see 
legend). The error bar in the upper-left corner represents the mean uncertainties of 14% in and 30% in S. The results of 

RDIO, a total of 329 clouds with M > 10^ Af 0 , are shown by the gray crosses. The mean results of |Solomon et aT] ( |l987| are also shown 
(open black square). The best-fit power-law relation oc is shown for both the IRDCs (n = 1.27 ± 0.52; dot-dashed line) and 

the GMCs {n = 0.70 ± 0 . 12 ; dashed line). Collectively, the IRDCs and GMCs have a best-fit power-law relation of n = 0.81 ± 0.10 (solid 
line). The dotted lines represents the expected scaling relation for virialized clouds with different values of avir = 0.5, 1, 2, 4, and 8 . 

uncertainty in g of ^ 14% based on 10% uncertainties in both a and R. 

Thus we considered the above scaling relationship for the results of our ten IRDCs and GMCs (Fig. |^. For IRDCs, 
S was estimated for each ellipsoidal cloud using the geometric mean observed radius R and the total ^^CO-derived 
mass. For GMCs, S was estimated using the circular area defined by the extraction radius R and its corresponding 
^^CO-derived mass. We estimate a random uncertainty in S of ^ 30% based on the uncertainty in Tex and potential 
variations in our assumed abundances. However, systematic uncertainties could be present at the factor of ^ 2 level. 

For the ten GMCs, we find that the clouds analyzed using gaussian profile fitting on the opacity corrected column 
density spectra to estimate the cloud mass and velocity dispersion, yield results most consistent with the scaling 
expected for gravitational virial equilibrium. Specifically, we find the best-fit power-law relation a jR}!"^ oc has 
n = 0.54 ± 0.05 and n = 0.70 ± 0.12, for the SE and CE extraction, respectively. We gauge the significance of these 
fits by estimating the Pearson’s correlation coefficient, p, and the one-sided p-value. Eor all the GMCs defined by CE 
and gaussian profile fitting on the opacity corrected column density spectra, we estimate p = 0.71 and a p-value of 
6.0 X 10“^, indicating a significant linearity in the data. 

On the other hand, for the ten IRDCs analyzed using gaussian profile fitting on the opacity corrected column 
density spectra to estimate the cloud mass and velocity dispersion, yield steeper best-fit power-law relations, with 
n = 1.57 ± 0.80 and n = 1.27 ± 0.52, for the SE and CE extraction, respectively. These power-law fits are less 
significant than those seen for the GMCs, with p = 0.56 and a p-value of 0.047 for the “CE,r,G” method. 

Eocusing on the “CE,r,G” extraction method, these IRDC and GMC results are displayed together in Eig. 
Overall the IRDCs and GMCs together yield a best-fit power-law relation of n = 0.81 ± 0.10, somewhat steeper than 
the n = 1/2 expected for virialized clouds. If a power-law relation with n = 0.5 is fit, then the normalization implies 
divir = 1-60 and 1.03 for IRDCs and GMCs, respectively, and dvir = 1-15 for the whole sample (with a Pearson’s 
correlation coefficient of p = 0.78 with a p-value of 1.7 x 10“^). We note that some of the dispersion in the data is 
due to cloud-to-cloud variation, which might be expected if there is systematic variation in degree of surface pressure 
confinement, magnetization, distance uncertainty (affecting R) or ^^CO abundance. 
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We conclude that the relatively tight correlation shown by the data in Fig. is evidence of cloud virialization 
over a wide range of scales within GMCs, perhaps representing a self-similar hierarchy of self-gravitating virialized 
structures. There is tentative evidence that IRDCs deviate from this virialized scaling relation in the sense of being 
more kinematically disturbed. However, the level of deviation of the IRDCs is quite modest and their best-fit dvir — 1.6 
is a value that is quite compatible with models of virialized singular polytropic {kp ^ 1.5) spherical clouds embedded 
in a pressure confining ambient medium. Such a model would yield dvir = 3.75 (but quite sensitive to kp). Including 
support from large-scale magnetic fields that lead to Alfven Mach numbers of unity would lower this to dvir = 1.34. 
However, given these uncertainties in kp^ together with those from abundance, using the absolute value of dvir 

to assess the degree of virialization is problematic, and we argue it is better to focus on the degree of correlation in 


the G versus S diagram. 

For comparison, the results of RDIO for a total of 329 clou ds with M > 10"^ 


2013 


13| 

Thi 


0 and with their S values scaled by 
are also shown in Fig. 


10 


The RDIO 


a factor of 0.49 to reflect our adopted abundance (see Tan et al. 

GMCs also follow the expected virialized scaling, but with larger scatter 
see larger scatter with our simpler cloud extraction methods 
have values of dvir — 1.1. 

We also note here that the absolute values of S that we derive for GMCs are similar to those of H 09 and RDIO 
H09 found an average mass surface density of S ^ 42 Mq pc“^ 


is is to be expected given that we also 
Consistent with our CMC results, the RDIO GMCs also 


( 1987|) ^^CO-defined 
2ior these larger 
> 10^ M0 yield 


using 162 of the Solomon et al. 

molecular clouds (but they noted that the expected values are likely underestimated by a factor 01" 
scale GMCs due to sub-thermal excitation of ^^CO). Additionally, the 329 RDIO clouds with M ^ lu iwiq 
S 88 Mq pc“^. For the 10 GMCs studied here, we find Sg — 80 pc“^ assuming CE, opacity corrected column 
densities, fitted gaussian line profile to derive mass, and circular cloud radii defined by For the 10 IRDCs, defined 
by the Simon et al. (2006) ellipses, we find Eg — 217 Mq pc“^. However, it is important to note that there is a real 
spread in E values between and among the considered IRDC and CMC scales of greater than a factor of 10. 


5. CONCLUSIONS 

We have presented a detailed study of the ^^CO gas on various size scales out to radii of 30 pc around 10 well-studied 
IRDCs (i.e., the sample of BT09 and BT12). We find that all the IRDCs are embedded in a ^^CO-defined CMC, i.e., a 
molecular cloud with M > 10^ Mq. We have taken care to assess the effects of cloud definition, specifically comparing 
simple extraction of a slab in position-velocity space with extraction only of connected emission features. We have 
also taken notice of the effects of multiple Gaussian components when they are present. We have considered column 
density estimates via assuming, first, optically thin ^^CO(l-O) emission (with excitation temperatures estimated from 
associated ^^CO emission) and then correcting for line optical depth effects. Our preferred method of “CE,r,G” cloud 
definition involves connected extraction, opacity-corrected column densities and fitting a single Gaussian profile to 
derive velocity dispersion and mass. 

We discussed the particular morphologies of the individual clouds in both position-velocity space and in relation to 
plane of sky imaging of the IRDCs. We find most IRDCs are offset from their 30-pc-scale CMC center by ^ 2 — 10 pc, 
but have only small velocity offsets of ^ 1km s“^. This suggests IRDCs have been formed by contraction/compression 
of structures that were already pre-existing with the GMCs rather than being formed from material swept-up by a 
shock. Six of the ten clouds have relatively simple morphologies in position-velocity space, while the other four exhibit 
more complex, multi-component structures. 

We then focussed on the kinematic and dynamical properties of the clouds. We measured velocity gradients as a 
function of size scale, 5, finding they decrease approximately as i.e., consistent with observations that IRDCs 

show larger velocity gradients than GMCs. Overall, this is also consistent with the gradients being caused by turbulence 
with a velocity dispersion versus size relation of the form a oc Such a scaling is seen in our sample of GMCs, but 
a somewhat shallower relation is apparent when considering IRDCs and GMCs together. Using velocity gradients to 
define projected rotation axes we do not see correlations with Galactic plane orientation and have approximately equal 
fractions of pro- and retrograde rotators with respect to Galactic orbital rotation (though this analysis is limited by 
the small sample of only 10 clouds). We measure rotational energy to gravitational binding energy ratios of P 0.045, 
with a slight increase seen on going from CMC to IRDC scales. 

Assessing virial equilibrium, we find cloud definition has a significant impact on the virial parameter, i.e., the ratio 
of internal kinetic to gravitational energy. CE,r,G-defined GMCs have dvir — 1.0, while IRDCs have dvir — 1.6. We 
argue that it is difficult to assess the dynamical state, i.e., how close to virial equilibrium, from the absolute values 
of dvir, especially given the uncertain effects of pressure confinement, large-scale magnetic field support and mass 
surface density estimation from ^^CO. However, these values dvir ^ 1.0 to 1.6 are quite compatible with Turbulent 
Core/Clump/Cloud models of singular polytropic spherical clouds that are pressure-confined by an ambient medium 
and have moderate large scale H-field support that yields Alfenic Mach numbers of order unity. Note, that ^^CO- 
defined clouds are likely to have significant surrounding molecular material that could provide pressure confinement. 

Rather than focus on absolute values of dvir, we follow H09 in examining the predicted scaling in the versus 

E plane. The results can be sensitive to the method of cloud definition. Eor our preferred method of CE,r,G-defined 
clouds, the GMCs, analyzed over a range of scales, show a strong correlation in this diagram, with a power law slope 
of oc E’^ with n = 0.70 ± 0.12, similar to the n = 1/2 value expected under virial equilibrium. Potentially this 

represents a self-similar hierarchy of self-gravitating virialized structures. 

However, the IRDCs exhibit a somewhat steeper scaling with n = 1.27 ± 0.52 (or n = 0.81 ± 0.10 for the combined 
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GMC and IRDC sample). This may indicate that IRDCs have moderately more disturbed kinematics than GMCs. 


et ah 1986 


which is qualitative exp ectation of models in which star formation is triggered by GMG-GMG collisions (e.g., Scoville 
' Tan||2QQQ ). However, it is also possible that systematic uncertainties affecting mass surface density 


Hernandez et ah 


2011 


2012), are playing a 


estimation from ^^GO, such as increased GO freeze-out in IRDGs (e.g., 
role. Improved methods of measuring mass and assessing degree of GO freeze-out in these clouds are needed to help 
distinguish these possibilities and thus reliably measure the dynamical state of gas that is likely to form star clusters. 

Additional future work will involve applying the techniques developed here to larger samples of clouds to better 
study connections with Galactic dynamics and Galactic environment, and potential systematic differences of GMGs 
that are forming IRDGs or that have relatively high dense gas mass fractions. 
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APPENDIX 

A. CLOUD MAPS: B-J 
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Fig. Al.— Same as Fig. but for Cloud B. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 20.7 — 31.6 km s“^. 
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Fig. A2. — Same as Fig. but for Cloud C. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 66.3 — 85.7 km s“^. 
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Fig. A3.— Same as Fig. but for Cloud D. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 84.1 — 93.7 kms“^. 
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Fig. A4. — Same as Fig. but for Cloud E. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 66.5 — 89.8 km s“^. 
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Fig. A5.— Same as Fig. but for Cloud F. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 54.3 — 65.6 km s“^. 
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Fig. A6.— Same as Fig. but for Cloud G. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 39.7 — 51.8 km s“^. 
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Fig. A7.— Same as Fig. but for Cloud H. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 39.7 — 47.1 km s“^. 
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Fig. A8.— Same as Fig. but for Cloud L The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 37.0 — 45.1 km s“^. 
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Fig. A9.— Same as Fig. but for Cloud J. The ^^CO(l-O) GRS integrated intensity map of the IRDC emission profile is over the 
velocity range: ulsr = 18-7 — 27.1 km s“^. 
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